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ABSTRACT 


tligh frequency (1IF) ionospheric propagation predictions have been available on 
mainframe computers since the late 1960s. Since the advent of low cost, computationally 
powerful personal computers, several propagation codes have been ported from 
mainframes to PC’s. This study compares results from two versions of the IONCAP and 
PROPHET HF propagation prediction codes to a database of measured electric field 
strengths. 

IONCAP-PC 2.5 predictions were compared to IONCAP-VAX 85.04 predictions 
and to those from PROPHET 3.2 (PC-version). A database of measured signal strengths 
from the CCIR containing over 16,000 points or cases was used as a benchmark for 
comparing code results. 

For both IONCAP-PC 2.5 and PROPHET 3.2, field strength predictions were low 
in more than 50°o of the cases examined. This was particularly true for PROPHET. 
which is considered to be a conservative model. PROPHET features quick solutions, 
graphical outputs, and a user-friendly environment. in comparison to IONCAP. 
JONCAP-PC 2.5, which is an improved version of a previous mainframe TONCAP, 
produced shghtiy more accurate predictions than IONCAP 85.04, and substantiails 


better results than those from PROPHET. 
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I. INTRODUCTION 


lonospheric propagation prediction programs are used daily in the Navy for HF 
communications planning. Thereforc, an assessment of the usefulness and accuracy of 
these programs is needed. In order to accomplish this, measured values must be com- 
pared with values calculated from the propagation prediction programs. Typically, three 
parameters can be compared: 

a. Maximum Usable Frequency (MUF), 

b. Signal-to-noise ratio (S N) at the receiver input terminals, 

c. Signal field strength at the receiver location. 

Measured and calculated MUFs have been compared in the past. Fortunately, the 
MUF is a preliminary parameter that is used for more complex calculations. 

Field strength comparisons between measured and calculated values are more pop- 
ular than S N comparisons, since estimating noise levels is not involved. 

In the real world during the propagation of a wave through the ionosphere, the 
following happens: There are numerous propagation paths possible between a transmit- 
ter and a receiver, depending on the vertical beamwidth of the transmitter’s antenna and 
the ionospheric lavers. Cven though a radio wave enters the ionosphere with a certain 
angle. it is split into two components (ordinary and extraordinary) because of the pres- 
ence of the earth’s magnetic field. Conditions in the ionosphere are not identical for the 
two components, Their reflections take place at different altitudes and their delay times, 
polarization and losses are different. Thev are also received at diflerent in-coming angles 
at the receiver. Therefore, more than one signal arrives at the receiver at the same time. 
Lach one has different amplitude, phase and polarization. Calculation of signal field 
strength at the receiver is a verv complicated process due to all of these factors. 

In order for the existing prediction programs to simulate the propagation through 
the ionosphere, the following are assumed: 

1. There is no interference between signals transmitted by the same source that ar- 
rive at the receiver through different paths; 

2. There is no birefringence (wave component splitting); 

3. The wave's polarization does not change by passing through the ionosphere, or 
after ground reflection and; 


4. There is no laver fluctuation during propagation. 





A. HISTORICAL BACKGROUND 

For many vears, numerous organizations have used the HIF spectrum to communi- 
cate over long distances. In the late thirties, communication svstems were subject to 
large Variations in performance. Effectiveness of long-distance systems increased in pro- 
portion with the ability to predict variations in the ionosphere, since such an ability 
pernutted the selection of optimum frequencies and other circuit parameters. To meas- 
ure 1onospheric parameters, a worldwide network of ionosondes was established. Vari- 
ations in signal amplitudes were observed and then recorded over various lif’ paths. 
Worldwide noise measurements were also taken. The results of this research showed that 
most of the variations in IIF path performance were directly related to changes in the 
ionosphere. which in turn were affected by solar activity, seasonal and diurnal variations. 
Latitude and longitude of the transmitter and the receiver also influenced performance. 

By 1948 a treatise of ionospheric radio propagation was published by the Central 
Radio Propagation Laboratory (CRPL) of the National Bureau of Standards, which 
outined the state-of-the-art in HF propagation. Manual techniques were given for an- 
alvzing HIP circuits of short, intermediate, and long distances. [Ref 1] 

Because these manual techniques were laborious and ume consumung. Various com- 
puter programs were developed to predict HIF circuit performance. All these programs 
were based on manual methods and used various numeric representauons of the 
ionospheric data. 

From 1945 to 1960, the Radio Propagation Laboratory (RPI) of the U.S. Army 
collected data concerning monthly median receiver input voltages for several circuits in 
the HF band. This data was statistically analvzed in order to study the parameters that 
influence propagation through the ionosphere, and a number of reports presented engi- 
neering methods of evaluating skv-wave field strengths [Ref. 2]. 

Recently, another bank of field strength data, called the CCIR HF data bank, was 
collected. This data bank contains 16,268 data points of observed electric field strength 
correlated with information for the month, year, hour (UT), frequency, and sunspot 
number and includes 181 combinations of circuit and frequency. The frequencies are 
from 2.5 to 26 MIlz, path Iengths are between 175 and 26,000 km, and the data spans 
the period of 1964 to 1985 for two sunspot cycle variations. 

The median monthly skvwave electric field strength (observed) is in one-hour time 
blocks. Its value is given in dB relative to ] wY m fora 1 KW ERP source. 


Professor A. Tomko from Johns Hopkins University used this data bank to make a 


statisucal analvsis of the prediction accuracy of HF propagation models (Ref. 3]. One 








of the models that he analyzed was IONCAP, mainframe version 85.04 (Ref. 4]. The re- 
sults of his analvsis will be used to compare with other models. In order to make an 
evaluation he compared each model's predictions of the electric field strength, £,. with 
empirical data, £,. from the CCIR IIF data bank, by computing the prediction error, 
e=E,— E,. 


B. ANALYSIS 

In the next chapters, data from the same CCIR data bank will be used, and the same 
methodology that was mentioned above will be followed to analyze the performance of 
IONCAP-PC 2.5 and PROPHET. 

First, the two versions of LONCAP will be compared and then IONCAP-PC 2.5 
with PROPHET. Finally, an evaluation of IONCAP-PC 2.5 and PROPHET will be 
conducted for paths, defined as polar. The appendix also contains a short analvsis of 
another HIF propagation model, called MINIPROP. 


Il. THEORETICAL BACKGROUND 


A. SKY WAVE PROPAGATION 
Skv wave propagation is defined as radio wave reflections in which ionospheric re- 
flections dominate over other mechanisms. 
1. The Ionosphere 
a. Structure 

The ionosphere consists of several electrically conducting layers abcve the 
Earth's surface which play a very important role in the propagation of radio waves. 
Radiation from the sun is the primary scurce of ionizavion of these lavers and controls 
their intensity. The ionization rate at various altitudes depends upon the intensity of the 
solar radiation and the ionization efficiency of the neutral atmospheric gases. Since the 
sun’s radiation is progressively absorbed in passing through the atmosphere, the residual 
ionizing ability depends upon the length of the atmospheric path, and consequently 
upon the solar zenith angle. The maximum ionization rate occurs when the solar zenith 
angle is zero, but geographic, diurnal and seasonal variations in ionization density also 
occur. 

The ionosphere is divided into three layers, designated ID, E, and F. respec- 
tively. These lavers may be subdivided. From the vic .point of HF propagation, the 
L-and F-lavers act mainly as radio wave reflectors and permit long range propagation 
between terrestrial terminals. The D-laver acts mainly as an absorber causing atten- 
uation in the lower HIF range. 

b.  D-layer 

The D-laver, at an altitude of 50-90 km above the earth, is the lowest layer. 
Its electron density exhibits large diurnal variations. The D-layer appears after sunrise 
and attains its maximum value shortly after local solar noon. It then gradually breaks 
down, and disappears after sunset. There is also a seasonal variation with maximum 
values during sumuner. Short waves penctrate the D-laver, but they are attenuated while 
passing through the layer. 

«. E-layer 

The E-laver, at 90-130 km, is subdivided into “normal’ and “sporadic” 

lavers. The former has its maximum electron density at [10 km altitude. It displavs a 


strong solar zenith angle dependence with maximum density near noon and a seasonal 








maximum in summer. It also displavs a solar cycle dependence with maximum densitv 
at solar sunspot maximum. The normal E-laver is important for daytime HF 
propagation at distances less than 2000 km. 

The E, is known as the sporadic E-laver because it arises only occasionally 
at most latitudes. It is located 120 km above the Earth and, when present, causes higher 
than normal critical frequencies. Its occurrence is strongly latitude dependent. At high 
Jautudes, it is essentially a night-time phenomenon; at low latitudes, a daytime phe- 
nomenon. 

d. F-layer 

The F-laver extends upwards from about 130 km and is divided into the F1- 
and F2-lavers. The Fl-laver exists only during davlight hours and is located 170 to 220 
km above the earth. The F2-laver is located 225 to 450 km above the earth and usually 
exhibits the highest ciectron density, which ranges from 10” electrons/ny in daytime, to 
about 3 x 10” electrons/n¥ at night. The F2-laver cannot be described by a simple model 
since it is strongly influenced by winds. diffusion and other dynamic effects. It is the 
principal reflecting laver used for long-distance communications. Its height and 
ionization density vary diurnally, seasonally and over the sunspot circle. At night, the 
Pl-laver merges with P2-laver at a height of about 300 km. 

2. Maximum Usable Frequency 

Predicting the performance of an TIF skyvwave path is a complex problem. 
However, there 1s much that can be learned from examining the dependencies of some 
of the relevant characteristics, of Which the maximum usable frequency is particularly 
important. 

The maximum usable frequency (MUF) for anv transmission distance is calcu- 
lated from the critical frequency multiplied by a factor that is a function of the trans- 
nussion distance (Figure 1 on page 6). In practice, the MUF is not a sharp limit since 
some propagation is possible on frequencies greater than the MUF. This happens since 
neither the ground nor the ionosphere are smooth reflectors. 

The critical frequencies of the E- and F2-lavers, known as fE and f,F2 respec- 
uvely, are the highest frequencies that can be reflected from the two layers at vertical 
incidence with the 1onosphere: they are related to the maximum electron densitics in 
those lavers. The value of (12 is usually greater than the value of fE because the 


electron concentrauion in the F-laver is ucuclly greater than in the [-laver. 
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Figure I. Range dependence of MUF factor for E- and F-layers. From [Ref. 5: p. 
64}. 


The plasma frequency fy is given by 
10.5 
I=IN, (2.1) 
where N is the electron density in electrons per cubic meter and fy 1s in Hertz. If N,, is 
the maximum clectron density in a given layer, then all waves for frequencies less than 


the plasma frequency that enter the ionosphere at vertical incidence will be reflected 


back to carth. The critical frequency is given by 
0.5 
fp = INO. (2.2) 
a. Variation of the critical frequency 


(1) Solar Cycle dependence. Figure 2 on page 8 shows the monthly mean 


noon critical frequencies of the E- and ['2-layers over the last solar cvele. It can be seen 


that there is a correlation of the F2-laver critical frequencies with the eleven-vear solar 
evcle. 

Additionally, the seasonal variation of £E is in phase with solar zenith 
angle (greater in the summer months). In contrast, the seasonal variation of f,F2 is in 
antiphase with solar zenith angle. This is known as the “winter anomaly”. Therefore, 
ualike the E- and the Fl-layer, the F2-laver does not follow a simple solar zenith angle 
either diurnally or seasonally. 

‘2; Annual, Seasonal, and Diurnal yariations. There is an annual 
variation of the monthly median (F2 . The greatest variation over the solar cvcle for 
davtime conditions occurs in winter and the least variation in the summer. By contrast, 
night-time conditions in the winter are mostly unaffected by sunspot number. 

There is a seasonal variation of the monthly median (F2. The largest 
Variation over the \ear occurs during periods of high sunspot activity. 

There is also a diurnal variation of the monthly median (F2. Diurnal 
Variations are most marked during winter months. 

3. Path Attenuation and Attainable Field Strengths 
For a specific shortwave path, when determining the transmitter power and an- 
tenna gains, the calculations are based on stable properties of the path. With propa- 
gauion conditions continuously varving. these path properties are characterized by 
attenuation and field strength values that can be expected for a given percentage of the 
time. 
a. Losses in a shortwave system 


If a shortwave link is considered as a closed transmission svstem, starting 


with the transmitter output and ending with the receiver input, its total losses L,,,, can 
be given from the following individual losses: 
Los = Lr + L,+ Lat Yeo G,- Gr (2.3) 


where 
L,, 1s the free-space attenuation in dB, 
L, is loss in the ionosphere in dB, 
L, is loss from ground reflections in dB, 
}; is a fading reserve in dB, and 
G, and G;, are the gains of the transmitting and the receiving antennas in 


JB, respectively. 
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Figure 2. Solar dependence of the critical frequency: © Monthly mean of the noon 
critical frequencies observed at Slough 1970-81 with superimposed 
sunspot activity. From [Ref. 5: p. 68]. 


(1) Free space attenuation. The decrease in radiated power density from 
a transmit antenna is a function of the square of the distance, with field strength de- ‘ 
creasing linearly with distance. Considering isotropic radiators for transmitting and re- 


ceiving, free space attenuation is given in dB from the relation: 


Ly, = 20 log(4ndgy] A) = 32.44 + 20 log £ + 20 log dey (2.4) 


where fis the operating frequency in MI Iz and d,, is the effective path length in km. 
(2) lonospheric losses. Radiowaves passing through the ionosphere suf- 
fer from loss due to absorption in the D-layer. This loss depends on the zenith angle 


of the sun, w, and the sunspot number R. 








The absorption index is given from the relation: 
1=(1 4+ 0.0037R){ cos 0.881)? (2. 


The absorption index is higher during winter months than during the rest of the vear, 
therefore, correction factors are used to account for the winter anomaly. 

The loss also depends on the elevation angle of the radiation and the 
effective frequency f+ fy. The effective frequency is the sum of the operating frequency 
and the gvrofrequency f,. The gvrofrequency is created by the Earth’s magnetic field 
(fy = 2.84 x 10" x B where B is the magnetic induction in Wb m’) . 

With absorption index, the effective frequency f+ f,,, and the angle 
of elevation, the loss in the ionosphere ZL, can be determined with the aid of the 
nomogram in igure 3 on page 10. 

‘3, Ground reflections loss. fa sky wave is reflected by the ionosphere 
back to earth, it reaches the ground with a polarization differing from that of the emitted 
wave. The energy of the incident wave will be split into vertical and horizontal compo- 
nents. Assuming that these two components are equal, the reflection loss from the 
ground, in dB. 


ROER, 
La= 10 log( —=—* ) (2.6) 


where R, is the reflection coefficient of the vertical component and R, is the reflection 
coefficient of the horizontal component. 

(d, Fading reserves. The losses described thus far refer to a stable radio 
path. However, field strength at the receiving point is influenced by continuously fluc- 
tuating propagation conditions along the path, due to several types of fading phenom- 
ena. The probability of occurrence of a specific field strength value follows a Rayleigh 
distribution. A fading reserve of 14 dB is, in general, recommended for shortwaves in 


order to attain availability 90% of the time. 





Losses per hop 
Frequency f +f, 





Example: Angle of elevation = 8”, absorption index f= 1.08. The line (1) connecting these (wo 
values intersects with the reference line at point P. The connecting tine (2) for the frequency 22.2 
Mlifz through P shows a foss of 6 dB per hop. 


Figure 3. Nomogram for determining the ionospheric losses per hop as a function 
of the absorption index, the elevation angle and the effective frequency. 


From [Ref. 6: p. 217]. : 








Ill. MODEL DESCRIPTIONS 


A. ITONCAP 

The lonospheric Communications Analysis and Prediction Program (IONCAP) is 
a program that attempts to solve the physics of each parameter [Ref. 7]. Therefore. its 
solutions are accurate but computationally intensive. For this reason, IONCAP is often 
used for communication circuit planning. 

IONCAP was designed in modular form, therefore, any subsection can be replaced 
Without affecung the rest of the program and is written in FORTRAN. The program 
is divided into seven independent subroutines: 

1. Input 

2. Path geometry 


. Antennas 


ft. Ue 


. lonospheric parameters 
5. Maximum usable frequency 


6. Svstem performance 


~I 


. Output. 

Input subroutines handle various input options. Three kinds of input are: card im- 
age, a long term data tape image, and an antenna tape image. Path geometry subrou- 
tines determune circuit geometry, select optimum areas to sample the ionosphere, and 
evaluate the magnetic field at these sample areas. Antenna subroutines process antenna 
data input cards and calculate antenna gains. 

lonospheric parameter subroutines evaluate all ionospheric parameters needed by 
the program. A complete explicit electron density profile is used in this program. The 
profile includes a D-E layer, an F2-laver, an Fl-ledge and an E-F vallev. 

Maximum Usable Frequency (MUF) subroutines evaluate MUFs and Optimum 
Traffic Frequencies (FOTs). The method used is a direct evaluation of the [requency that 
Propagates 50%o of the time for the MUF and 90°%o of the time for the FOT. The E-, 
Fl- and F2-laver MUFs are computed. The sporadic-E MUF is also computes. 

System performance subroutines evaluate all the needed circuit parameters. There 
are two different sets of subroutines. one for shorter distances (less than 10,000 km) and 


one for longer distances. However, the user may run cither set of subroutines for anv 


given distance (overriding the 10,000 km breakpoint). 








The short path model evaluates each possible ray path for the circuit; high and low 


angle, E, Fl, and F2 modes, over the MUF modes and sporadic-E modes. Losses include 
regular D-E absorption and sporadic losses. For lower frequencies which have low re- 
flection heights (less than 90 km), a correction to the frequency dependence parameter 
was added. 

Bv using the antenna gains and the ionosphere-loss function, the model evaluates 
the optimum radiation angle at the transmitter and the optimum reception angle at the 
receiver. The 29 output subroutines generate all output options as line images. 

The IONCAP program performs four basic analysis tasks. These tasks are summa- 
rized below: 

1. lonospheric Parameters 

The ionosphere is predicted by using parameters which describe the four 
ionospheric lavers: E, Fl, F2, and E,. At each sample area the location, time of day, and 
all ionospheric parameters are derived. These can be used to find an electron density 
profile. By integration, a plot of reflection height against frequency can be constructed, 
which is called an ionogram. 

2. Antenna Patterns 

The user may precalculate the antenna pattern and gain used for system per- 
formance predictions. 

3. Maximum Usable Frequency (MUF) 

The maximum frequency at which a skvwave mode exists can be predicted. The 
90% (FOT), 50° (MUF), and 10% (Ilighest Probable Frequency) levels are calculated 
for each of the four ionospheric lavers. These numbers are a description of the state of 
the 1onosphere between two locations and not a statement on the actual performance 
of any operational communication circuit. 

4. Svstem Performance 

A comprehensive prediction of radio system performance parameters (up to 22) 


is provided. 


B. PROPHET 

The Advanced PROPHET HF prediction is the follow-on successor to the 
SOLRAD PROPHET family of real-time HF propagation prediction and assessment 
systems (Ref. 8 and 9]. It has been designed to give conservative and real-time solutions, 


in order to support US military tactical comniunications. 


In order to achieve this. PROPHET does not attempt to solve the physics of each 
parameter, like most of the classical IIF prediction models do. PROPHET solves the 
problem bv using simplified empirical models which have been thoroughly verified with 
HIF oblique incidence data (ionograms). This provides rapid soluuions with small and less 
expensive computers. 

The Advanced PROPHET system consists of over 24 ionospheric algorithms and 
several special application models. 

The basic model set used in Advanced PROPHET is contained in a module called 
“PROFLIB”. This module calculates the various parameters, like the MUF, and then 
applies them to an application routine which produces the final product. 

The first quantity needed for a given skywave channel is the channel bandwidth or 
frequency spread between the MUF and the Lowest Usable Frequency (LUF). An al- 
gorithm known as MUF-S8S has been developed in order to calculate the MUP. A 
comparable algorithm to calculate the LUF also has been developed. The LUI model 
is complicated by absorption effects. 

On several of the PROPHET products, an additional frequency, the Frequency of 
Optimum Transmission (FOT). is noted. The FOT is defined as the frequency with 90" 
reliability, and is given as a guide to the user in selecting the optimum frequency for the 
current mussion. , 

The bandwidth for a given channel depends primarily on solar effects in the 
ionosphere. Hlowever. communicating information over this channel depends not only 
on the frequency, but also on various enviromental effects such as noise, absorption, 
required signal-to-noise ratio for the class of emission and antennas. .A model for the 
expected noise for various environments has been developed. It includes contributions 
from atmospheric, galactic and man-made sources. An algorithm calculates the 
ionospheric absorption index for any time and location. System parameters and antenna 
patterns are input to PROPHET by the user. An algorithm converts geographic coor- 
dinates to range between transmitter and receiver and azimuth beam headings. 

An algorithm computes the ray paths that will propagate between the designated 
transmitter and receiver, for the specified operating frequency and for a specific model 
of the ionosphere. The calculations are controlled by the initial launch angle, launch 
angle increment, and maximum launch angle used in computing the ray paths. The cal- 
culations start with the initial angle and compute successive ravs for angles increased 
by an increment until two rays are fot nd that bracket the receiver site. An interpolation 


routine is then used to find the launch angle that will hit the receiver. The resulting ray 
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path is called a mode. The procedure continues until the launch angle reaches the 
maximum linut or the number of specified modes is reached. Influences by a jammer site 
are included, and models are used to assess the effects they have on the receiving svstem. 

Short-range propagation bv ground or surface wave depends on several factors such 
as frequency, antenna height, ground conductivity, dielectric constant. and polarization. 
Empirical models are used to simulate this propagation. 

The field strength model was derived by comparing the output of the model with an 
extensive set of measurements which depends on the MUF, LUF, the effective radiating 
power and environmental parameters. It provides a method for rapid assessement of 
HIF skvwave signal quality as a function of frequency within the propagation bandwidth 
for links between specified points on the Earth’s surface. It also provides matrices of the 
signal-to-(noise + jammer) ratio as functions of time and frequency. 


The above models comprise the basis for HF modeling. as expressed in PROPHET. 


I4 








IV. CODIPARISON OF IONCAP 85.04 WITH IONCAP-PC 2.5 


A. ASSUMPTIONS 
Statistical analysis of the prediction accuracy of IONCAP 85.04 has been conducted 
by Professor A. Tomko and the results are contained in Reference 3. Comparison of 
IONCAP 85.04 predictions with CCIR HF data bank measurements encompassed all 
values (16,268) documented in the CCIR data bank. Values from TIONCAP-PC 2.5, 
which is an improved version of IONCAP 85.04, were compared to a limited number of 
values from the CCIR data bank. 
Three factors were considered that mainly influence the propagation through the 
ionosphere: 
1. The Path Leneth 
A short path was considered to be a path less than 2500 km. A medium path is 
a path greater than 2500 km but less than 8000 km. A long path is a path greater than 
SQ0U km. The same classification of paths has been done in Reference 3. Therefore, 
short. medium, and long paths were chosen for IONCAP-PC 2.5 in the same proportion 
to the paths that were contained in the CCIR data bank. 
2. The Season of the Year 
The data were classified in two periods, winter and summer. Winter was con- 
sidered October to March and summer from April to September for the northern hemi- 
sphere. For a fair sampling, the data was chosen in the same proportion as in the CCIR 
for these two periods. 
3. The Sunspot Number 
The sunspot number was an important factor for choosing the sample data. The 
CCIR data was sorted by category: long, medium, and short paths. Each category was 
then plotted versus the sunspot number. The samples were chosen so that their distrib- 
ution versus sunspot number was similar to the distribution of the CCIR data but ona 
smaller scale. 
The lines of the graphs, obtained from the analysis of IONCAP-PC 2.5, merely 


connect data points. 





B. ANALYSIS 
1. Macroscopic View of All Path Data 

Figures 4 and 5 show the frequency distribution of the predicted errors for the 
two versions of JONCAP. The standard deviation of the error distribution of 
IONCAP-PC 2.5 is a bit smaller than the standard deviation of IONCAP 85.04. How- 
ever, this does not imply better performance for IONCAP-PC 2.5 because the median 
of the error distribution for IONCAP 85.04 is closer to zero than the median error of 
IONCAP-PC 2.5. 

Figures 6 and 7 show the cumulative distribution of the absolute value of pre- 
diction errors for the two versions of IONCAP. Comparing Figure 6 with Figure 7, it 
is seen that 8S°o of LONCAP-PC 2.5 data is included between -20 dB and + 20 dB error, 
while for IONCAP 85.04, only 80°o of the data is included. Therefore. it is concluded 
that IONCAP-PC 2.5, in general, is slightly more accurate than IONCAP $5.04 for this 
measured data base. 

In Figures § and 9, the frequency distribution of the prediction errors for the 
two Versions of IONCAP is shown for various path lengths. Figures 10 and 11 show the 
cumulative distribution of the absolute value of the prediction errors for the two versions 
of IONCAP, for various path lengths. The peak values of the error frequency for 
JONCAP-PC 2.5 are higher (35%, 48%, and 31% for short, medium, and long paths 
respectively) than the respective peak values of TONCAP 85.04 (13%, 12.5% 0 and 6.3%). 
The short path data represents 4$°%o of all the data. The error value that corresponds 
to the peak value of error frequency for short path data is 7 dB for IONCAP-PC 2.5 and 
2 dB for [ONCAP 85.04. For long and medium path data, the error values that corre- 
spond to peak values of error frequency have no significant difference for the two ver- 
sions of TONCAP. Therefore, because of the significant difference in error value of short 
path data, it is noted that. while there is such a significant difference in peak values of 
error frequency, there is not a significant difference in the final performance of the two 
different versions of IONCAP. 

2. Paths Comparison 

IONCAP 85.04 predictions are more accurate for short and medium paths than 
for long paths as shown in Figure 10. Figure 11 shows that for IONCAP-PC 2.5, pred- 
ictions are a bit more accurate for short and medium paths than for long paths. For 
IONCAP 85.04 (Figure 10), it can be concluded that the random component of predic- 


tion error increases with increasing path length. The same cannot be said for 
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Figure 4. Frequency distribution of [ONCAP 85.04 prediction errors. From (Ref. 
3). 
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Frequency distribution of LONCAP-PC 2.5 prediction errors. 
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Figure 6. Cumulative distribution of the absolute value of LONCAP 85.04 prediction 
errors. From [Ref. 3]. 
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Figure 7. Cumulative distribution of the absolute value of LONCAP-PC 2.5 predic- 


tion errors. 
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Figure 8. Frequency distribution of IONCAP 85.04 prediction errors for various 
. path lengths. From [Ref. 3}. 
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Figure 9. Frequency distribution of FONCAP-PC 2.5 prediction errors for various 


path lengths. 
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Figure 10. Cumulative distribution of the absolute value of IONCAP 85.04 predic- 
tion errors for various path lengths. From [Ref. 3]. 
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Figure 11. Cumulative distribution of the absolute value of TONCAP-PC 2.5 pre- i 





diction errors for various path lengths. 








JONCAP-PC 2.5 since the prediction error is less for medium paths than for short paths 
(Figure 11). 


Figures 12, 13, and 14 show in detail, the statistics of the short, medium and 


long paths, respectively, for IONCAP-PC 2.5. These figures also reflect that field 
strength predictions for IONCAP-PC 2.5 are consistently low on medium and long 
paths, since the median error for both paths is negative (-6 and -8 dB respectively). The 


same conclusion is contained in Reference 3 for IONCAP 85.04. 
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path lengths, for 1.102 samples. 
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Figure 13, Frequency distribution of LONCAP-PC 2.5 prediction errors for medium 


path lengths. 


27 


ee 
= 
m 
O 
mY 
rm 
1) 


= -8 dB 
= -10.5 dB 
= 17.7 dB 


652 


Median Error 


Samples 
Mean Error 
Std Dev Error 


vay S WM ro) 
foe) ice) rt 


Quasjeg) RONANOAYA YOUN 





Tigure 14. Frequency distribution of LONCAP-PC 2.5 prediction errors for long 
path lengths. 


V. COMPARISON OF IONCAP-PC 2.5 WITH ADVANCED PROPHET 


A. ASSUMPTIONS 

The same CCIR data bank, used in Reference 3, was used to make a statistical 
analysis for prediction accuracy of IONCA” ind PROPHET. For making a comparison 
of these two programs, 1.961 cases were run using Method 17 for IONCAP and the Field 
Strength option for PROPHET. The lines of the graphs contained in this chapter merely 


connect data points. 


B. ANALYSIS 

Medium and long path data for IONCAP-PC 2.5, that has alreadv been analyzed in 
Figures 13 and 14, respectively, will be used for the comparison with PROPHET. The 
results of the analvsis for LONCAP-PC 2.5, that are presented in Figure 15 will be used 
for the comparison with PROPHET, for short path data. 

1. Short Path Data 

The error frequency versus error for PROPHET for the short paths is shown in 

Figure 16. It is seen in Figures 15 and 16 that the median error of IONCAP is closer 
to zero than the error from PROPHET. The standard deviation error for IONCAP is 
much smaller than the error from PROPHET. Therefore, IONCAP is much more ac- 
curate than PROPHET for short path predictions. The same is concluded from Figure 
17, which shows the cumulative distribuuon of the absolute value of prediction errors. 
for short paths, for both models. Note that 90% of IONCAP data is included between 
-20 dB and +20 dB error, however for PROPHET, only 65% of the data is included in 
this range. 

2. Medium Path Data 

In Figure 18, the error frequency versus the error is shown, for PROPHET, for 

medium paths. From Figures 13 and I$ it is difficult to conclude which of the models 
is more accurate, since their median error values and standard deviations are close. 
However, it can be seen from Figure 19 that 92°0 of IONCAP data is included between 
-20 and + 20 dB error, while for PROPHET only 80%o is included. Therefore, IONCAP 
is more accurate than PROPHET for medium paths. 

3. Long Path Data 


In Figure 20, the error frequency versus the error is shown, for PROPHET, for 


long paths. In Figure 21, the cumulative frequency versus error is shown for both pro- 
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Figure 15. Frequency distribution of LONCAP-PC 2.5 prediction errors for short 


paths, for 783 samples. 
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Figure 16. Frequency distribution of PROPHET prediction errors for short paths. 
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Figure 17, Cumulative distribution of the absolute value of TONCAP and 
PROPHET prediction errors for short paths. 
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Figure 18. Frequency distribution of PROPITET prediction errors for medium 


paths. 
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Figure 19. Curnulative distribution of the absolute value of TONCAP and 
PROPHET prediction errors for medium paths. 
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grams, for long paths. By comparing Figure 14 with Figure 20 it is seen that the median 
error 1s very close to zero for PROPHET (-1 dB), and relatively high for IONCAP (-8 
dB). The standard deviation is a little higher for PROPHET (19.9 dB) than for IONCAP 
(17.7 dB). Therefore, it can be concluded that PROPHET is slightly more accurate than 


IONCAP for long path predictions. This is more obvious in Figure 21 where 82% of 
PROPHET data is included between -20 dB and +20 dB of error while only 78% of 
1ONCAAP data is included. 

4. Macroscopic View of All Path Data 

The prediction analysis for IONCAP and PROPHET for the data from all paths 
is presented in Figures 22 and 23 respectively. It can be seen that the two models have 
about the same median error. The standard deviation error is much higher for 
PROPHET. Therefore, IONCAP is more accurate than PROPHET. In Figure 24, the 
cumulative frequency versus the absolute error is shown for IONCAP and PROPHET 
for data from all paths. It can be seen that 81% of IONCAP data is included between 
-20 dB and +20 dB error and only 70°o of PROPIIET data is included. In general, it 
is concluded that IONCAP is more accurate than PROPHET. 

5. Paths Comparison for PROPHET 

In Figure 25, error frequency versus error is shown for short, medium and long 
paths. for PROPHET. This information was accumulated from Figures 16, 18 and 20. 
In Figure 26, cumulative frequency versus absolute error is shown, for these paths, for 
PROPHET. 

It is seen in Figures 25 and 26 that PROPHET yields high errors during pred- 
ictions of short path circuits. 

By comparing the median error of PROPIIET, for short, medium and long paths 
(Figures 16, 18. and 20), it is concluded that PROPHET field strength predictions are 
consistently low on short and long paths, as could be expected for a conservative pro- 
gram. 

6. Polar Paths 

For this study, a path is defined as polar if, 

a. One Or both stations of a circuit are located above 60° latitude, or 

b. any part of the path between stations crosses the zone between the 60° and 
90° latitude. 

Almost all of the polar paths of this CCIR data bank have both stations below 
60° latitude, and a portion of the path crosses the zone between the 60° and 90° latitude. 
Unfortunately, there is not any path in the CCIR data bank that crosses the Auroral 
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Figure 20. Frequency distribution of PROPHET prediction errors for long paths. 
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Figure 21. Cumulative distribution of the absolute value of IONCAP and 
PROPHET prediction errors for long paths. 





Figure 22. 
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Frequency distribution of IONCAP prediction errors. 
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Figure 23, 





Frequency distribution of PROPIIET prediction errors. 
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Figure 24. Cumulative distribution of the absolute value of IONCAP and 
PROPHET prediction errors. 
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Frequency distribution of PROPHET prediction errors, for short, me- 


dium and long paths. 
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Figure 26. 
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Cumulative distribution of the absolute value of PROPHET prediction 
errors, for short, medium and long paths. 


zones, Where severe ionospheric disturbance occurs. The “polar” paths that are analyzed 
below are all classified as long paths and their stations are below 60° latitude. 

The frequency distribution of IONCAP and PROPHET prediction errors is 
shown in Figures 27 and 28 respectively, for polar paths. It is known that most 
ionospheric models are less accurate for predictions of polar paths than for other paths. 
For PROPHET, this can be observed, since the standard deviation for the polar paths 
is the highest of those already seen (23.4 dB). However, this is not true for JONCAP, 
since the median error of polar paths (-5 dB) is closer to zero than the median error of 
long paths (-8 dB) which already contain polar paths. The standard deviation of polar 
paths (15.3 dB) is also less than the standard deviation of long paths (17.7 dB). There- 
fore, IONCAP is better at predicting the field strength of polar paths, as previously de- 
fined, than predicting long paths, as supported by Figures {4 and 27. Figures 11 and 29 
also demonstrate that, for IONCAP, between -20 and +20 dB error, 86% of the data 
of polar paths is contained. On the other hand, only 76°%o of the data of long paths is 
contained. This also supports the fact that the paths, defined as polar, are not actually 
auroral region polar paths. 

In Figure 29, the cumulative distribution of the absolute value of the prediction 
error is shown, for IONCAP and PROPHIET, respectively. It can be seen that 86% of 
IONCAP data is included between -20 and + 20 dB error while only 73°¢ of PROPHET 
data is included. Therefore, IONCAP is more accurate than PROPHET for “polar” path 


predictions. 
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Figure 27. Frequency distribution of IONCAP prediction errors for polar paths. 
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Figure 28. Frequency distribution of PROPHET prediction errors for polar paths. 


45 


< 
O 
Z 
iS) 
| 


--- PROPHET 





i=) 
o oO jon) S&S So i] =) 
On ive) é \O -y vt ioe) 


100 


QuadJeg) AONANOAUA AAILWINWNOD 


Figure 29. Cumulative distribution of the absolute value of LONCAP and 


PROPHET prediction errors, for polar paths. 
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ABSOLUTE VALUE OF ERROR (dB) 


VI. CONCLUSIONS AND RECOMMENDATIONS 


. The conclusions that were extracted from the last two chapters are summarized be- 
low when comparing to the CCIR field strength data base: 
? e IONCAP 85.04 predictions are most accurate for short, and medium paths [Ref. 
3h. 
¢ IONCAP-PC 2.5 predictions are most accurate for medium paths. 


e The random component of the prediction error for IONCAP 85.04 increases with 
increasing path length [Ref. 3]. This is not true for the prediction error of 
IONCAP-PC 2.5. 


e IONCAP 85.04 [Ref. 3] and IONCAP-PC 2.5 field strength predictions are con- 
sistently low on medium and long paths. 


@ TONCAP-PC 2.5 field strength predictions are, in general. slightly more accurate 
than IONCAP 85.04 field strength predictions. 


e Jn IONCAP-PC 2.5, a great improvement is observed in predicting long path cir- 
cults In comparison with JONCAP 85.04. The improvement for predictions of 
medium paths circuits is small, and there is no improvement for predictions of short 
path circuits. 


e The random component of the prediction error for PROPHET increases with de- 
creasing path length. since PROPHET is less accurate for predictions of field 
strength of short path circuits than IONCAP. 


; e JONCAP is generally more accurate than PROPHET. 
¢ PROPHET ficld strength predictions are consistently low on short and long paths. 
i * TONCAP is more accurate than PROPHET for “polar” path field strength pred- 
ictions. 

In the future. a further analysis of MINIPROP should be conducted. A model 
called AMBCOM has become available at NPS [Ref. 10]. The special features of this 
code are a model of the electron density profile in the high latitude ionosphere (including 
the auroral zone) and a model for computing the auroral absorption. PROPHET 4.0 
is available and a new version of IONCAP-PC for polar regions (ICECAP) will also be 
available soon. An evaluation of these models must also be conducted. 

An area of great interest concerning 1onospheric propagation, is the area near the 
poles. Professor Tudor Jones from the University of Leichester in England has collected 
data for measured electric field strength for true polar paths (Ref 11]. An analvsis should 


be conducted by using this data bank for ©1] models. 
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APPENDIX MINIPROP RESULTS 


MINIPROP 3.0 is a PC-based computer program for predicting ionospheric propa- 
gation which covers the HF band (3-30 MHz). It predicts received signal levels using an 
algorithm that searches the potential ionospheric modes to find the mode that provides 
the strongest received signal level [Ref. 12]. 

IONCAP and PROPHET outputs are electric field strength predictions in dB above 
1 uV.m, while the primary MINIPROP output is a voltage prediction at the input of the 
receiver in dB above 0.5 uV. Therefore, MINIPROP output must be transformed to field 
strength values in order to be compared with the measured values of the CCIR data 


bank by using the formulas: 


F? 


a" 320z 





where 
P, is the power flux at the antenna in Watts m~, 


E is the electric field strength in Volts,m; and 


where 0 
P is the received power at the input of the receiver in Watts, 
V is the voltage at the input of the receiver in Volts, 
R is the receiver input resistance in Qhms (R = 5022), 
A, is the effective aperture of the antenna in m’. 


For a halfwave dipole antenna the effective aperture A, is given by 
2 
2 
A,= Pais 1.64 


where / is the wavelength in meters. 


By using the formulas mentioned above, the electric ficld strength E is given by 


E=V+20logf—32 
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where 

E is the electric field strength in dB above I nV/m, 

V is the voltage at the input of the receiver in dB above 0.5 nV and 

fis the frequency in MHz. 

In order to evaluate the performance of MINIPROP, 122 cases were run with 
IONCAP, PROPHET and MINIPROP. All these cases were for long paths. 

In Figure 30, the frequency distribution of the prediction errors for MINIPROP is 
shown. In Figure 31, the cumulative distribution of the prediction errors for the three 
models is shown. From Figure 30 it is seen that the median error for MINIPROP is 14 
dB. For the same data, the median errors for IONCAP and PROPHET are -9 dB and 
-10 dB, respectively. The standard deviation error for MINIPROP is 15.2 dB; for 
IONCAP and PROPHET 15 dB and 20 dB, respectively. It is concluded that while 
IONCAP and PROPHET give conservative predictions, MINIPROP gives optimistic 
predictions. From Figure 31, it is concluded that MINIPROP has about the same per- 
formance as PROPHET and is less accurate than IONCAP for the cases analyzed. 


However, since only one category of paths was analvzed and the number of the 


samples was limited, a further study for the performance of MINIPROP is needed. 
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Figure 30. Frequency distribution of MINIPROP prediction errors. 
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Figure 31. Cumulative distribution of the absolute value of LONCAP, PROPHET 
and MINIPROP prediction errors. 
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